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ABSTRACT 

We explore the possibility of setting stringent constraints to the Dark Energy equation 
of state using alternative cosmic tracers like: (a) the Hubble relation using HII galaxies, 
which can be observed at much higher redshifts (z^> 3.5) than those currently traced by 
SNIa samples, and (b) the large-scale structure using the clustering of X-ray selected 
AGN, which have a redshift distribution peaking at z ~ 1. 

In this paper we use extensive Monte-Carlo simulations to define the optimal strat- 
egy for the recovery of the dark-energy equation of state using the high redshift (ztt 2) 
Hubble relation, but accounting also for the effects of gravitational lensing, which for 
such high redshifts can significantly affect the derived cosmological constraints. We 
investigate the size of the sample of high-z HH-galaxies needed to provide useful con- 
straints in the Dark Energy equation of state. Based on a "Figure of Merit" analysis, 
we provide estimates for the number of 2^ z'tS 3.5 tracers needed to reduce the cosmo- 
logical solution space, presently provided by the Constitution SNIa set, by a desired 
factor. The analysis is given for any level of rms distance modulus uncertainty and we 
find that an expected reduction (i.e. by ~ 20% — 40%) of the current level of HH-galaxy 
based distance modulus uncertainty does not provide a significant improvement in the 
derived cosmological constraints. It is much more efficient to increase the number of 
tracers than to reduce their individual uncertainties. 

Finally, we propose a framework to put constraints on the dark energy equation 
of state by using the joint likelihood of the X-ray AGN clustering and of the Hubble 
relation cosmological analyses. A preliminary joint analysis using the X-ray AGN 
clustering of the 2XMM survey and the Hubble relation of the Constitution SNIa set 
provide: fi m = 0.31±0.01 and w= -1.06±0.05. We also find that the joint SNIa-2XMM 
analysis provides significantly more stringent cosmological constraints, increasing the 
Figure of Merit by a factor ~ 2, with respect to that of the joint SNIa-BAO analysis. 



1 INTRODUCTION 



We live in a very exciting period for our understanding of the 
Cosmos. Over the past decade the accumulation and detailed 
analyses of high quality cosmological data (eg., supernovae 
type la, CMB temperature fluctuations, galaxy clustering, 
high z clusters of galaxies, etc.) have strongly suggested that 
we live in a flat and accelerating universe, which contains at 
least some sort of cold dark matter to explain the clustering 
of extragalactic sources, and an extra component which acts 
as having a negative pressure, as for example the energy of 
the vacuum (or in a more general setting the so called dark 
energy), to explain the observed accelerated cosmic expan- 



sion (eg. Riess, et al. 1998; 2004; 2007, Perlmutter et al. 
1999; Spergel et al. 2003, 2007, Tonry et al. 2003; Schuecker 
et al. 2003; Tegmark et al. 2004; Seljak et al. 2004; Allen et 
al. 2004; Basilakos & Plionis 2005; 2006; 2009; Blake et al. 
2007; Wood-Vasey et al. 2007, Davis et al. 2007; Kowalski et 
al. 2008, Komatsu et al. 2009; Hicken et al. 2009, Amanullah 
et al. 2010, etc). 

Due to the absence of a well-motivated fundamental 
theory, there have been many theoretical speculations re- 
garding the nature of the dark energy (hereafter DE), on 
whether it is a cosmological constant or a field that provides 
a time varying equation of state, usually parametrized by: 
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PQ = w(z)pq , 



(i) 



with pq and pQ the pressure and density of the exotic dark 
energy fluid. For a large class of DE models, we have: 



w(z) = wo + wi/(z) , 



(2) 



with wo =w(z = 0) and f(z) an increasing function of red- 
shift. A particular example of w(z) is its 1st order Taylor's 
expansion around w(0), which provides f(z) — z/(l + z), ie., 
the so-called CPL form of the DE equation of state (Cheva- 
lier & Polarski 2001, Linder 2003; see also Peebles & Ra- 
tra 2003, Dicus & Repko 2004; Wang & Mukherjee 2006). 
Of course, it could also be conceived that the equation of 
state parameter does not evolve cosmologically (quintessence 
Dark Energy model; QDE). 

It is clear that one of the most important questions 
in Cosmology and cosmic structure formation is related 
to the nature of dark energy (as well as whether it is 
the sole explanation of the observed accelerated expan- 
sion of the Universe) and its interpretation within a fun- 
damental physical theory (eg., Albrecht et al. 2006; Pea- 
cock et al. 2006). To this end, a large number of very 
expensive experiments are proposed and/or are at var- 
ious stages of development, viz the Dark Energy Sur- 
vey: http://www.darkenergysurvey.org/, the Joint Dark 
Energy Mission: http://jdem.gsfc.nasa.gov/, HETDEX 
http : //www . as . utexas . edu/hetdex/, 

Pan-STARRS: http://pan-starrs.ifa.hawaii.edu, Eu- 
clid: http://sci.esa.int/euclid/, Wfirst: 
http : //wfirst . gsf c .nasa. gov, etc. 

Therefore, the paramount importance of the detection 
and quantification of DE for our understanding of the cos- 
mos and for fundamental theories implies that the results 
of the different experiments should not only be scrutinized, 
but alternative, even higher-risk, methods to measure DE 
should be developed and applied as well. It is within this 
paradigm that our current work falls. Indeed, we wish to 
constrain the DE equation of state using, individually and 
in combination, the Hubble relation and large-scale struc- 
ture (clustering) methods, but utilizing alternative cosmic 
tracers for both of these components. 

Thus, we will trace the Hubble relation using HII galax- 
ies, which can be observed at higher redshifts than reached 
by current SNIa surveys to distances where the Hubble re- 
lation is more sensitive to the cosmological parameters. The 
HII galaxies can be used as standard candles (Melnick, Ter- 
levich & Terlevich 2000, Melnick 2003; Siegel et al. 2005; 
Plionis et al. 2009) due to the correlation between their ve- 
locity dispersion, Hp luminosity and metallicity (Melnick 
1978, Terlevich & Melnick 1981, Melnick, Terlevich & Moles 
1988). Furthermore, the use of such alternative high z tracer 
will enable us to check the SNIa results and lift any doubts 
that arise from the fact that they are the only tracers of the 
Hubble relation used to-date (for possible usage of GRBs 
see for example, Ghirlanda et al. 2006, 2009; Basilakos & 
Perivolaropoulos 2008 and references therein) *. 



* GRBs appear to be anything but standard candles, having a 
very wide range of isotropic equivalent luminosities and energy 
outputs. Nevertheless, correlations between various properties of 
the prompt emission and in some cases also the afterglow emis- 
sion have been used to determine their distances. A serious prob- 



Additionally, we use X-ray selected AGN at a median 
redshift of ~ 1, which is roughly the peak of their redshift 
distribution (see Basilakos et al. 2004; 2005, Miyaji et al. 
2007), in order to determine their clustering pattern and 
compare it with that predicted by different cosmological 
models (see Matsubara 2004). 

Although each of the previously discussed components 
of our project (Hubble relation using HII galaxies and angu- 
lar/spatial clustering of X-ray AGN) will provide interest- 
ing and relatively stringent constraints on the cosmological 
parameters, especially under our anticipation that we will 
reduce significantly the corresponding random and system- 
atic errors, it is the combined likelihood of these two type 
of analyses that enables us to break the known degeneracies 
between cosmological parameters and determine with great 
accuracy the DE equation of state (see Basilakos & Plionis 
2005; 2006; 2009). 

Below we present the basic methodology and expecta- 
tions of the two components of our method. In section 2 we 
present the details of the first component where we develop 
a Monte-Carlo simulation approach designed to ultimately 
provide a rule of thumb of how many HII galaxies we need to 
obtain a particular level of the DE equation of state param- 
eter uncertainty. We also develop a method to account for 
the effects of gravitational lensing, which at such high red- 
shifts are significant. In section 3 we present the details of 
the second component and in section 4 we present an exam- 
ple of joining the two components to provide cosmological 
constraints. The conclusions are listed in section 5. 



2 COSMOLOGICAL PARAMETERS FROM 
THE HUBBLE RELATION 

In the matter dominated epoch and in flat universes, the 
Hubble relation depends on the cosmological parameters via 
the following equation: 



H(z) = H E(z) 
with 

e 2 (z) = n m (i+ z f+n k {i+z) 2 +n Q exp 



Jo 1 + x 



(3) 



dx 



which is simply derived from Friedman's equation. We re- 
mind the reader that fi m , fl k and Qq(= 1 — fim — are 
the present fractional contributions to the total cosmic mass- 
energy density of the matter, the spatial curvature and dark 
energy source terms, respectively. 

Supernovae SNIa are considered standard candles at 
peak luminosity and therefore they have been used not only 



lem that hampers a straightforward use of GRBs as Cosmological 
probes is the intrinsic faintness of the nearby events, a fact which 
introduces a bias towards low (or high) values of GRB observ- 
ables and therefore the extrapolation of their correlations to low- 
2 events is faced with serious problems. One might also expect 
a significant evolution of the intrinsic properties of GRBs with 
redshift (also between intermediate and high redshifts) which can 
be hard to disentangle from cosmological effects. Finally, even if a 
reliable scaling relation can be identified and used, the scatter in 
the resulting luminosity and thus distance modulus is still fairly 
large. 
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Figure 1. Left Panel: The expected distance modulus difference between the DE models shown and the reference A-model (w= — 1) 
with f2 m = 0.27. Right Panel: The expected distance modulus differences once the f2 m -w(z) degeneracy is broken (imposing the same 
Q m value as in the comparison model). 



to determine the Hubble constant (at relatively low red- 
shifts) but also to trace the curvature of the Hubble relation 
at high redshifts (see Riess et al. 1998, 2004, 2007; Perlmut- 
ter et al. 1998, 1999; Tonry et al. 2003; Astier et al. 2006; 
Wood-Vasey et al. 2007; Davis et al. 2007; Kowalski et al 
2008; Hicken et al. 2009; Amanullah et al. 2010; Wang, Li 
& Li 2011; Kim 2011; March et al. 2011; Adak, Bandopad- 
hyay, Majumdar 2011). In practice one relates the distance 
modulus of the SNIa to its luminosity distance, dr,, through 
which the cosmological parameters enter: 



H = m- M = 5 log d L + 25 



and 



dh = sin 1 1 



i?0V"fc 



da; 
o W) 



(5) 



(6) 



which for a flat universe (ilk ~ 0) reduces to: 



c(l + z) 
H 



da; 



(7) 



The main result of numerous studies using this procedure is 
that distant SNIa's are on average dimmer by ~0.2 mag than 
expected in an Einstein-deSitter model, which translates in 
them being ~ 10% further away than expected. 

The amazing consequence of these results is that we live 
in an accelerating phase of the expansion of the Universe, an 
assertion that needs to be scrutinized on all possible levels, 
one of which is to verify the accelerated expansion of the 
Universe using alternative to SNIa's extragalactic standard 
candles. Furthermore, the cause and rate of the acceleration 
is of paramount importance, ie., the DE equation of state is 
the next fundamental item to search for and to these direc- 
tions we hope to contribute with our current project. 



2.1 Theoretical Expectations: 

To appreciate the magnitude of the Hubble relation vari- 
ations due to the different DE equation of states, we plot 
in Figure 1 the relative deviations of the distance modulus, 
An, of different dark-energy models from a nominal standard 
(w= —1) A-cosmology (with Q m = 0.27 and Qa ~ 0.73), 
with the relative deviations defined as: 



Afl = (t*A ~ Mmod 



(8) 



The parameters of the different models used are shown in 
Figure 1. As far as the dark-energy equation of state pa- 
rameter is concerned, we present the deviations from the 
standard model of two models with a constant w value and 
of two models with an evolving equation of state parame- 
ter, utilizing the form of eq.(2). In the left panel of Figure 
1 we present results for selected values of £l m , while in the 
right panel we use the same dark-energy equations of state 
parameters but for the same value of Q m (= 0.27) (ie., we 
avoid the Q m — w(z) degeneracy). 

Three important observations should be made from Fig- 
ure 1: 

(i) The relative magnitude deviations between the dif- 
ferent dark- energy models are quite small (typically ^ 0.1 
mag), which puts severe pressure on the necessary photo- 
metric accuracy of the relevant observations. 

(ii) The largest relative deviations of the distance moduli 
occur at redshifts z£ 1.5, quite larger than those currently 
traced by SN la samples, and 

(iii) There are strong degeneracies between the different 
cosmological models at redshifts 1, and in some models 
even up to much higher redshifts (eg., between the mod- 
els with (fi m ,w ,wi) = (0.31,-1,0) and (0.29,-1,0.3); see 
Figure 1). 

Luckily, as discussed already in the introduction, such degen- 
eracies can be broken by using other cosmological probes (eg. 
the clustering of extragalactic sources, the CMB shift pa- 
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Figure 2. Probability density function of the lensing magnification (left), flux (middle) and magnitude (right) distributions for a source 
with intrinsic magnitude uncertainty of a- lnt = 0.1 at two different redshifts (thin line corresponds to z = 1, while the thick line to z = 3). 



rameter, BAO's, etc). Indeed, current evidence overwhelm- 
ingly show that the total matter content of the universe is 
within the range: 0.2J; f2 m i$ 0.3, a fact that reduces signif- 
icantly the degeneracies between the cosmological parame- 
ters. 



2.2 Gravitational Lensing Effects on High z 
Distance Moduli 

A potentially important systematic effect that could hin- 
der attempts to put stringent cosmological constraints via 
the Hubble relation, especially when using high z standard 
candles (which as we saw are those precisely that differen- 
tiate between DE equations of state) is related to gravita- 
tional lensing by structures intervening between source and 
receptor. It is indeed known that the gravitational poten- 
tial of large-scale structure affects the propagation of light 
from high redshift sources and thus also the distance mod- 
ulus of similarly high redshift standard candles (eg., Holz & 
Wald 1998; Holz & Linder 2005; Brouzakis & Tetradis 2008 
and references therein). These studies assume a Robertson- 
Walker background superimposing a locally inhomogeneous 
universe and take into account both strong and weak lensing 
effects. The resulting magnification distribution of a single 
source over different paths is non- Gaussian and therefore 
has a non-trivial effect on its distance modulus, especially 
so for the high z standard candles. 

The main characteristics of the magnification proba- 
bility density function, P(fi a ), as derived from a variety 
of studies based on Monte-Carlo analyses and ray-tracing 
techniques, is that P(jJ. a ) resembles a log- normal distribu- 
tion with zero mean (the mean flux of each source over all 
possible different paths is conserved, since lensing does not 
affect photon numbers), with a mode shifted towards the 
de-magnified regime and a long tail to high magnification. 
This implies that most sources will be de-magnified, induc- 
ing an apparently enhanced accelerated expansion, while a 
few will be highly magnified. The effect is obviously stronger 
for higher redshift sources since the lower the redshift the 
less the optical depth of lensing. Note that although the de- 
tailed shape of P(/J. a ) is a function of the underlying cosmol- 
ogy, density profile and evolutionary phase of the intervening 



cosmic structures, the main features discussed previously re- 
main unaltered (eg., Wang et al. 2002). 

We will therefore model the lensing effect using a log- 
normal magnification distribution, according to appendix A 
of Holz & Linder(2005; HL05 hereafter). The fact that the 
mean flux, over all different paths of a source, converges to 
the unlensed value implies that if we had a large number 
of standard candles densely populating all the redshift bins, 
the lensing effects would be smoothed out and it would be 
unnecessary to correct. However, this is not usually the case 
and therefore we need to take lensing into account (espe- 
cially for the high z sources). 

Two main effects of lensing will be accounted for: 

• the increase of the distance modulus uncertainty by a 
further term due to lensing, <T e ff, which was found by HL05 
to be a linear function of redshift with a e g = 0.093z. For a 
large number, N, of paths (or equivalently of sources) the 
lensing distribution is approximately Gaussian with variance 
a% and although the lensing distribution of a single path 
(source) is non- Gaussian, we can define the effective variance 
of a single path (or source) as: = Na% . As suggested by 
HL05 a reasonable /N contribution to the total distance 
modulus variance, is given by requiring N'Si 10 within z-bins 
of ~ O.lz width. Note that this is the only lensing dependent 
effect that has been taken into account in some of the SNIa 
based analyses (eg., Kowalski et al. 2008; Amanullah et al. 
2010). 

• the shift of the mode of the distance modulus distribu- 
tion to de-magnified values (fainter) due to lensing. This is 
an effect that has not yet been taken into account in the 
SNIa based studies. 

In order to investigate this later effect, and as we have al- 
ready pointed out previously, we will use the log-normal ap- 
proximation to the magnification distribution due to its nice 
analytical properties and its resemblance to the actual mag- 
nification distribution (see HL05). If [t, a is the source mag- 
nification, then its probability distribution is approximated 
by: 



P<J*a) 



2n SfJ, a 



exp 



(In fig - (In^g)) 2 
2S 2 



(9) 
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Figure 3. The mode offset (<5m) and effective variance (cr e g) (a) of a source as a function of redshift (left panel) and (b) as a function 
of the number of sources, TV, tracing the same redshift (for this example z = 3). 



with S the variance of In/ia- The mean magnification is 
given by (/i a ) = exp ((In fj, a ) + S 2 /2) = 1, implying that 



(In fj,a) 



-S /2. Therefore, the probability function is 



skewed ((In jj, a ) < 0) and determined by only the S parame- 
ter. From eq.(9) HL05 derived the corresponding flux distri- 
bution, which is also log-normal, and then the corresponding 
magnitude distribution, which is given by: 



P{m) = 



exp 



(— m — m + b(ln^ a )) 2 
2^1 



(10) 



with b = 2.5/ In 10 and a 2 n = a 2 bs + (bS) 2 . Therefore, we 
re-capitulate that the effects of lensing are: 
(a) an offset of the mean, given by: 



(m) = mo + 5m 



(11) 



where mo is the intrinsic (de-magnified) magnitude, and 
8m = — b(ln/j,a} = bS 2 /2, and, 

(b) an increase of the variance for which we have that a 2 s = 
{bS) 2 and thus S 2 = a 2 a /b 2 . 

Recalling that for large enough sources (N), in relatively 
small z-bins (~ O.lz), we have that a 2 s = Na 2 ^, we obtain 
that the magnitude offset of sources within the redshift bin 
is given by: 



_ _ (0.093,) 2 
v ' 2bN 2bN 
and the total source variance: 



o2.00 = ^obs + (0.093z) 2 /N 



(12) 



(13) 



In Figure 2 we plot as an example the magnification, 
flux and magnitude probability distributions for the case of 
a single source placed at two different redshifts [z = 1 and 
3) and for an intrinsic observational magnitude uncertainty 
of o"obs = 0.1. The intrinsic flux of the source is assumed to 
be 10 -17 erg s _1 cm -2 and the intrinsic apparent magnitude 



m = 18. The main effects discussed previously are clearly 
seen, ie., a mode offset towards the demagnificaton (fainter) 
regime and an enhanced variance, both increasing with red- 
shift. In Figure 3 we plot the expected increase of both the 
mode offset and variance (in magnitudes) as a function of 
redshift for a single source (left panel). In the right panel 
of the same Figure we plot the suppression of both quan- 
tities as we increase the source sampling (the case shown 
corresponds to a source located at a redshift 2 = 3). 

We will therefore correct statistically the distance mod- 
uli of observed standard candles (SNIa, GRBs, HII-galaxies, 
etc) by subtracting an offset 8m(z) from their raw distance 
modulus (according to eq.12), within redshift bins of ~ Q.lz 
width and using as the total distance modulus uncertainty 
that given by eq.(13). 



2.3 Best Strategy to Determine the DE Equation 
of State 

2.3.1 Fitting Models to the Data 

We can now proceed with our investigation to find an effi- 
cient strategy to put more stringent constraints on the dark- 
energy equation of state. To this end we have decided to 
re-analyse two recently compiled SNIa samples, the Davis 
et al. (2007) [hereafter D0T\ compilation of 192 SNIa (based 
on data from Wood-Vasey et al. 2007, Riess et al. 2007 and 
Astier et al. 2007) and the Constitution compilation of 397 
SNIa (Hicken et al. 2009). Note that the two samples are 
not independent since most of the D07 is included in the 
Constitution sample. 

Firstly, we present in the left panel of Figure 4 the Con- 
stitution SNIa distance moduli overploted (red-continuous 
line) with the theoretical expectation of a flat cosmology 
with (fi m ,w) = (0.27,-1). In the inset we plot the dis- 
tance moduli difference between the SNIa data and the pre- 
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Figure 4. SNIa distance moduli as a function of redshift. Inset 
Panel: Distance moduli difference between the best fit model (see 
Table 1) and the SNIa data. The blue dashed line is the corre- 
sponding difference between the reference (red-continuous line) 
(Q m ,w)= (0.3, -1.04) and the (« TO ,w)= (0.3, -0.85) dark-energy 
models. 
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Figure 5. Cosmological parameter solution space using either of 
the two SNIa data sets (Constitution: red shaded contours and 
DOT. black contours). Contours corresponding to the 1 and 3<r 
confidence levels are shown (ie., plotted where — 21n£/£ max is 
equal to 2.30 and 11.83, respectively). Inset Panel: Normalized 
redshift distributions of the two SNIa data sets (the shaded his- 
togram corresponds to the Constitution set). 



viously mentioned model. To appreciate the level of accu- 
racy needed in order to put constraints on the equation 
of state parameter, we also plot the distance moduli dif- 
ference between the reference (f2 m ,w) = (0.27, —1) and the 
(f2 m ,w) = (0.27, —0.85) models (continuous red and dashed 
blue line respectively). 

We proceed to analyse the SNIa data by denning the 
usual likelihood estimator^ as: 



£ SNIa (p) oc exp[-xsNia(p)/2] 



(14) 



where p is a vector containing the cosmological parameters 
that we want to fit for, and 



XSNIa(p) = ^ 



P)-M° bs (^) 



(15) 



where fj, th is given by eq.(5), Zi is the observed redshift and 
Gi is the distance modulus uncertainty, which includes the 
observational uncertainty and the gravitational lensing vari- 
ance (see eq. 13). Since in occasions the observational dis- 
tance modulus uncertainty has the form: [ij'j , ie., it is non- 
symmetric (due to its logarithmic dependence on the flux), 
we will use a slightly different weighting scheme in the min- 
imization function that takes into account the asymmetric 
observational uncertainty. Following Barlow (2004), and as- 
suming that the likelihood function of the observed n, de- 
rived from the theoretical /i(p), is a Gaussian, we can use 



the following weighting scheme of the \ function: 

r th / \ obsl 
<Ji — <7\ + (7 2 [/i (Pj — A* J 



t Likelihoods are normalized to their maximum values. 



(16) 



with <j\ = 2a p a n /(a p + a n ) and a?, = (a p — a n )/(<Tp + a n ). 
Obviously, when a p = a„ we recover the usual symmetric 
error weighting. 

In what follows we will constrain our analysis within 
the framework of a flat (Q m + = 1) cosmology and 
therefore p = (O m ,wo,wi). Note that we sample the var- 
ious parameters on a grid as follows: the matter density 
fi m G [0.04,0.64], the equation of state parameter w G 
[—2.0, —0.5], while when using a time-dependent equation 
of state: wo G [0,-2] and Wi G [—3,3]. The typical step 
size that we use is 0.0015. Note that the uncertainty of each 
fitted parameter will be estimated after marginalizing one 
parameter over the other, providing as its uncertainty the 
range for which A^ 2 < 2.3 (2cr). Such a definition, however, 
may hide the extent of a possible degeneracy between the 
two fitted parameters and therefore it is important to vi- 
sualize the 2D solution space, as indicated in the relevant 
contour plots. 



2.3.2 Larger Numbers? 

The first issue that we wish to address is how better have 
we done in imposing cosmological constraints by increasing 
the available SNIa sample from 181 to 366 (excluding the 

2 < 0.02 SNIa)^, ie., more than doubling the sample. Table 
1 presents various solutions using each of the two previously 



t We use only SNIa with z > 0.02 in order to avoid redshift 
uncertainties due to the local bulk flow (eg. Hudson et al. 1999 
and references therein). 
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Figure 6. Comparison between the Constitution SNIa con- 
straints (red shaded contours) and those derived by a Monte- 
Carlo procedure designed to closely reproduce them (for clarity 
we show only contours corresponding to 1 and 3 a confidence lev- 
els). Inset Panel: The Constitution SNIa distance modulus devi- 
ations from the best fit model (Q m , w) ~ (0.30, —1.01); see Table 
1 - and a random realization of the model deviations (red shaded 
histogram). 



mentioned samples. Note that since only the relative dis- 
tances of the SNIa are accurate and not their absolute local 
calibration, we always marginalize with respect to the in- 
ternally derived Hubble constant (for methods that do not 
need to a prion marginalize over the internally estimated 
Hubble constant, see for example Wei 2008). 

Regarding the fitted parameters uncertainty, we remind 
the reader that the definition we use (see above) cannot 
clearly reveal the extent of the degeneracy between the two 
parameters. A possible measure of such a degeneracy, be- 
yond inspecting the relevant contour plots, is to also present 
the whole range of the ler contours for each parameter. For 
example, the corresponding ranges are: Q m S [0.11,0.42] 
and Q m G [0.18,0.40] for the D07 and Constitution SNIa 
data sets respectively, while the corresponding w ranges are: 
wG [-0.66,-1.48] and wG [-0.78,-1.38]. 

Although the derived cosmological parameters are con- 
sistent between the two data sets, possibly indicating the 
robustness of the method, the corresponding goodness of fit 
(the reduced \ 2 ) ls significantly larger in the case of the 
Constitution set (1.21 compared to 1.045 of the D07 set). 
This appears to be the outcome of the different approaches 
chosen in order to join the different contributing SNIa sub- 
sets. According to Hicken 2009 (private communication) in 
the case of the D07 the nearby SNIa were constrained to 
provide a x 2 /dof ~ 1 by hand, while no such fine-tuning 
was imposed on the UNION set (on which the Constitution 
set is based). A secondary reason could be that the latter 
set includes distant SNIa which have typically larger dis- 
tance modulus uncertainties, with respect to those used in 
D07. Overall, the higher x 2 /dof value of the Constitution set 



should be attributed to a typically lower uncertainty in fi. 
As a crude test, we have increased by 20% the distance mod- 
ulus uncertainty of the Constitution nearby SNIa (z^ 0.4) 
and indeed we obtain x 2 /dof ~ 1.07, similar to that of D07. 
To also test whether lensing could have a significant effect on 
the derived cosmological values, we apply our lensing magni- 
fication correction procedure (see section 2.2) to both SNIa 
compilations and find a very small and insignificant change 
of the uncorrected for lensing results (see Table 1, last two 
rows), but interestingly a slightly better reduced \ 2 value. 

In Figure 5 we can also see that although the SNIa sam- 
ple has doubled in size, the well-known banana shape region 
of the (fi m ,w) solution space, indicating the degeneracy be- 
tween the two cosmological parameters, is roughly the same 
for both data sets. However, there is a reduction of the size 
of the solution space when using the Constitution SNIa com- 
pilation (see also Table 1) at roughly the level expected from 
Poisson statistics. 

A first conclusion is therefore that the increase by 
~ 100% of the Constitution sample has not broken the de- 
generacy in the (f2 m , w) plane and thus has not provided 
significantly more stringent constraints to the cosmological 
parameters. We have further verified that the larger number 
of SNIa's in the Constitution sample are not preferentially 
located at low-z's (see inset panel of Fig. 5) - in which case we 
would not have expected more stringent cosmological con- 
straints using the latter SNIa sample, but they have very 
similar ^-distributions. 

We already have a strong hint, from the previously pre- 
sented comparison between the D07 and Constitution re- 
sults, that increasing the number of Hubble relation tracers, 
covering the same redshift range and with the current level 
of uncertainties, as in the available SNIa samples, does not 
appear to be an efficient avenue for providing stringent con- 
straints of the cosmological parameters. 

2.3.3 Lower uncertainties or higher z's: 

We now resort to a Monte-Carlo procedure to investigate 
which of the following two directions, that bracket many 
different possibilities, provide the required more stringent 
cosmological constraints: 

• Reduce significantly the distance modulus uncertainties 
of SNIa, tracing however the same redshift range as the cur- 
rently available samples, or 

• use tracers of the Hubble relation located at redshifts 
where the models show their largest relative differences 
(z^ 2), with distance modulus uncertainties comparable to 
that of the highest redshift SNIa's ((<r M ) ~ 0.4). At such 
large redshifts however we expect that the gravitational lens- 
ing magnification/de-magnification effects will be significant 
and therefore we will also use the algorithm presented in sec- 
tion 2.2 to statistically degrade the intrinsic source flux and 
investigate its effects on the derived cosmological parame- 
ters. 

The Monte-Carlo procedure is based on using the ob- 
served high z SNIa distance modulus uncertainty distribu- 
tion (ct m ) and a model to assign random pi-deviations from 
a reference H(z) function, that reproduces exactly the orig- 
inal banana-shaped contours of the (fl m ,w) solution space 
of Figure 5, or in the case of the CPL model of the DE 
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Table 1. Cosmological parameter fits using the SNIa data for a flat prior cosmology. Note that for the case where p = (f2 m ,w), the 
indicated uncertainties are estimated by fixing one parameter at its best value and allowing the other to vary, providing as its uncertainty 
the range for which A% 2 < 2.3. 





D07 




Constitution 




w 




xL n /dof w 




xL n /dof 


Raw 


-1 

-1.005 ±0.076 


0.287 ±0.020 
0.289 ± 0.030 


186.721/180 -1 
186.721/179 -1.038 ±0.053 


0.285™ 
0.300 ± 0.022 


439.745/365 
439.703/364 


Lensing corrected 


-1 

-0.995 ± 0.075 


0.288 ± 0.020 
0.286 ± 0.030 


184.775/180 -1 
184.775/179 -1.036 ± 0.053 


0.284 ±0.014 
0.299 ± 0.022 


438.263/365 
438.229/364 
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Figure 7. Comparison of the model Constitution SNIa constraints (black contours) with those (filled contours) derived by reducing to 
half their uncertainties (left panel), with those derived by adding a sample of 76 high z tracers (2^5 z^S 3.5) with a distance modulus 
mean uncertainty of (<r M ) ~ 0.5 and no lensing degradation (central panel), and with those by including statistically the expected lensing 
degradation (right panel). For clarity we show only contours corresponding to the 1 and 3 a confidence levels. 



equation of state the corresponding contours in the u>o,Wi 
solution-space. Indeed, after a trial and error procedure we 
have found that by assigning to each SNIa (using their true 
redshift) a distance modulus deviation (Sfi) from a refer- 
ence model having a Gaussian distribution with zero mean 
and variance given by the observed (cr^) 2 , and using as the 
relevant indi vidual distance modulus uncertainty the follow- 
ing: of = U (1.2Sfii) 2 + 4> 2 (with cf> a random Poisson devi- 
ate within [—0.01, 0.01]) we reproduce exactly the banana- 
shaped solution range of the reference model. This can be 
seen clearly in Figure 6, where we plot the original Consti- 
tution SNIa solution space (red shaded contours) and the 
model solution space (black contours). In the inset panel 
we show the distribution of the true SNIa deviations from 
the best fitted model as well as a random realization of the 
corresponding model deviations. 

Armed with the above procedure we can now address 
the questions posed previously. Firstly, we reduce to half 
the random deviations of the SNIa distance moduli from 
the reference model (with the corresponding reduction of 
the relevant uncertainty, Oi). The results of the likelihood 



analysis can be seen in the left panel of Figure 7. There is a 
reduction of the range of the solution space, but indeed quite 
a small one. Secondly, we add to the Constitution SNIa sam- 
ple, a mock subsample of 76 high z tracers with a distance 
modulus mean uncertainty of (a^) ~ 0.5 (corresponding to 
that of the current HH-galaxy data) randomly distributed 
between 2^ z^ 3.5, ie., in a range where the largest devia- 
tions between the different cosmological models occur (see 
Figure 1) . We now find a significantly reduced solution space 
(central panel of Figure 7), which shows that indeed by in- 
creasing the H(z) tracers by a few tens, at those redshifts 
where the largest deviations between models occur, can have 
a significant impact on the recovered cosmological parameter 
solution space. If we include the expected lensing degrada- 
tion of the distance modulus (according to eq.12), then we 
observe (right panel of Figure 7) a slightly worsening of the 
solution space, but still significantly smaller than that of the 
left panel of Fig. 7. 

The main conclusion of the previous analysis is that a 
more efficient strategy to decrease the uncertainties of the 
cosmological parameters, based on the Hubble relation, is 
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Figure 8. The "reduction" parameter S, indicating the factor by 
which we reduce the 2a contour area of the cosmological param- 
eters (Q m ,w) solution space (QDE model) as a function of the 
number of high z tracers (2 < z < 3.5) of the Hubble relation 
and for two different values of the mean intrinsic distance modu- 
lus scatter (as indicated in the plot). Circular points correspond 
to using the high z tracers together with the current best SNIa 
data set, while the squares to using only the high z tracers (and 
a local z < 0.2 calibration sample). Inset Panel: The "reduction" 
parameter for the case of using 100 high z tracers as a function 
of the mean distance modulus uncertainty, (<Ta). The lines corre- 
spond to logarithmic fits to the data (see text). 



to use standard candles which trace also the redshift range 
2^$ zj$ 3.5. However, in such a case the effects of gravita- 
tional lensing can be severe, especially for small number of 
high z tracers, and therefore it is necessary to be taken into 
account . 



2.3.4 Figure of Merit Analysis 

In order to study the relation between the number of high 
z tracers used and the corresponding reduction of the cos- 
mological parameter solution space, we will use the Figure 
of Merit (FoM; Bassett 2005; Albrecht et al. 2006; Bueno 
Sanchez, Nesseris & Perivolaropoulos 2009), defined as the 
reciprocal area of the 2a contour (ie., where — 2bxC/£ max = 
6.14) in the parameter space of any two degenerate cosmo- 
logical parameters [eg., (fi m , w) for the QDE model or (wo, 
Wi) for the CPL model], and which has been found to be 
a useful measure of the effectiveness of a data set in con- 
straining cosmological parameters. A larger figure of merit 
indicates a greater accuracy in constraining the cosmological 
parameters. 

Here we will use a slightly different quantity, which we 
call "reduction factor" and is indicated by S, defined as the 
ratio of the FoM of the SNIa+high z Hubble relation solution 
to that of only the SNIa (in both cases we use the Consti- 
tution set), in order to study the question of how better can 
we constrain the cosmological parameter space, when adding 
iVhigh z high z tracers of the Hubble relation, with respect to 



the best current SNIa data set as a function of the number 
of high z tracers. For example a value 5 = 2 indicates that 
the FoM based on the SNIa+high- z Hubble relation is half 
of that based on the Constitution SNIa data set, ie., the 2a 
range of the solution space is reduced by a factor of 2. 

In Figure 8 we present the results of our analysis as 
a function of the number of high-z tracers for the QDE 
model. We present results only for the realistic lensing degra- 
dation case and for two different values of the distance 
modulus mean observational uncertainty, ranging between 
a pessimistic ((<7 M ) = 0.5; open circles) and an optimistic 
((<?>) = 0.25; filled points) value. 

It is evident that including even a small number of high- 
z tracers we can reduce significantly the cosmological pa- 
rameter solution space. There is a roughly linear relation 
between S and the number of high-z tracers used, JVhi g h z , 
which depends obviously on the distance modulus mean un- 
certainty, (<j m ) and on whether one uses in addition to the 
high-z tracers also the lower redshift SNIa data. 

It is also interesting to note that the high-z HII galaxies 
could constrain cosmological parameters with the level of ac- 
curacy provided by current SNIa data sets (for JVhi g h z ~ 200) 
and thus lift any doubts that arise from the fact that SNIa 
are the only reliable tracers of the Hubble relation used to- 
date. Of course the above relatively large number of high-z 
HII galaxies can be significantly reduced by including an 
intermediate population of HII galaxies, ie., tracing similar 
depths as the current SNIa samples (zi$ 1). In such a case 
the expected values of S will be intermediate between the 
only high-z and SNIa-high-z curves of Figure 8. For exam- 
ple, a realistic case of a sample with 80 0.2 < z < 1 and 60 
zzl 2 HII galaxies and a rms distance modulus uncertainty 
of ~ 0.35 will provide similar constraints as the current SNIa 
based analyses (S ~ 1). 

In order to quantify the previous results and provide 
a tool to estimate the number of high-z tracers necessary 
to reduce the current SNIa solution space by a given fac- 
tor, taking into account the whole parameter space, we first 
normalize the S values by that given for, say, A?high z = 100 
(Sioo). We then quantify how S100 depends on (cr M ), which is 
shown in the inset panel of Figure 8. The continuous curves 
are logarithmic fits to the data, which are given by the fol- 
lowing equations: 



qQDE 



1.871og 10 (( ( 7 M > 
0.441og 10 ((^> 



-1 
-1 



0.74) + 1.28 
0.15) +0.30 



SNIa+high z 
only high z 



(17) 



Then, in order to obtain the number of high z tracers, -/Vhi g h z , 
necessary to reduce the cosmological solution space (in the 
QDE model) by a factor S, we fit the normalized value, 
S/ S100 and find: 



iVhi 



highz 



187 SySioo 
139 SySioo 



39 



SNIa+high z 
only high z 



(18) 



which has a typical uncertainty of <7jv — ±5 for both the 
SNIa+high z and only high z cases. The continuous (red) lines 
in the left panel of Fig.8 are derived from eq.(17) and eq.(18) 
for (it,,) = 0.25 and 0.5, and it is evident that they repro- 
duce extremely well the observed S values (points). As an 
example, we can ask how many high z tracers, with say 
(cr M ) — 0.4, do we need to add to the current SNIa data 
set in order to reduce by a factor of 2 (S = 2) the current 
SNIa QDE solution space. Using eq.(17) and eq.(18) we find 
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A7high z — 80, which drops to ~ 60 for (cr M ) = 0.25. It is there- 
fore interesting to point-out that a reduction by a factor of 
2 in the distance modulus uncertainty, of the high z tracers 
(which is really a non-trivial aim) can be compensated by a 
relatively small increase in the number of high-z tracers. 

Repeating the previous analysis for the case of an evolv- 
ing DE equation of state (CPL; as in eq. 2), and after 
marginalizing with respect to Q m , we also find a reduction 
of the (wo, wi) solution space, when we include the high- 
z tracer subsample (Fig. 9), but significantly smaller than 
that of the QDE parametrization. We can again estimate 
what is the necessary number of high-z tracers, JVhi g h z , hav- 
ing a mean distance modulus error of (cr M ), in order to re- 
duce the cosmological (wo, wi) solution space by a factor S. 
Again using a parametrization based on the value of S for 
Miigh z = 100, we have that: 

cpl / 0.491og 10 «cr, 1 )- 1 + 0.65) + 1.09 SNIa+high, 



c 1 - 
'5 100 

and 



0.20 log 10 ((<7„) -1 + 0.69) + 0.22 only high z 



high z 



404 S/Sioo 
211 S/Sioo 



300 
106 



SNIa+high z 
only high z 



(20) 



with a typical uncertainty of <tjv — ±17 and ±7 for the 
SNIa+high— z and only high-z cases, respectively. These re- 
sults imply than in order to reduce by a factor of 2 (S = 2) 
the current SNIa CPL solution space using high-z tracers 
with (ct m ) = 0.35, one needs Nhigh z — 300 and 1200 for the 
SNIa+high-z and only high-z tracers case, respectively. The 
latter value can be significantly reduced if we include an 
intermediate redshift (0.2,$ z^S 1) HII sample, as discussed 
also for the QDE case previously. For example, using a sam- 
ple of 80 such intermediate z HII galaxies reduces this num- 
ber by a factor of 2. In any case, the large number of high-z 
tracers of the Hubble expansion, needed to effectively con- 
strain the CPL equation of state, renders this task rather 
unrealistic. Therefore, in order to provide stringent cosmo- 
logical constraints for the CPL model (ie., the values of wo 
and wi), it would be necessary (a) to combine the high- 
z Hubble relation with that of current SNIa data, and (b) 
to join the Hubble relation analysis with other cosmological 
tests, like the one that is an integral part of our proposal, ie., 
the clustering of X-ray AGN. Of course other cosmological 
probes, like BAOs, can and should be used as well. 

We now draw the main conclusions of our Monte-Carlo 
analysis: 



Even a small number of high-z (2 



3.5) tracers 



of the Hubble expansion can reduce significantly the QDE 
model parameter solution space. 

• For the case of the CPL model, in order to reduce the 
(wo, wi) solution space, provided by the current Constitu- 
tion SNIa set, by the same amount as in the corresponding 
QDM model, one needs three or more as many high-z Hub- 
ble expansion tracers. 

• It appears that the effort to reduce significantly the cur- 
rent level of random distance modulus scatter of HII galaxies 
is not as important as it is to increase the number of high-z 
HII galaxies, unless one is able to reduce it to (cr M )$ 0.1 — 0.2 
(as can be seen in the inset panel of Fig. 8). 



Ifi 



(19) 




200 



Figure 9. As in Figure 8 but for an evolving DE equation of state 
(CPL model) and after marginalizing with respect to Q m - The 
input cosmological model has (f2 m ,wo, wi) = (0.3, —0.98, —0.48). 
The axes scale has been kept as in Figure 8 in order to appreciate 
the significant reduction of efficiency of the high z Hubble relation 
in providing cosmological constraints for the CPL model. 



2.4 A High z Hubble Relation tracer: HII galaxies 

HII galaxies, compact extragalactic objects experiencing 
massive bursts of star formation, have a high luminosity 
per unit mass, in large part concentrated in a few strong 
emission lines in the optical rest frame. This ensures that 
the first, obvious, requirement for a standard candle to be 
usable at very large distances is met. 

The potential use of HII galaxies as distance indicators 
stems from the fact that as one increases the mass of the 
young stellar component, not only the ionizing output in- 
creases, but also the turbulent velocity of the gas, which 
is indicative of supersonic motions in the gas in the stellar 
gravitational potential, becomes larger. This effect induces 
a correlation between the integrated luminosity in a nebular 
hydrogen recombination line, e.g. L(H/3), which is propor- 
tional to the number of ionizing photons, and the line width 

(T. 

Terlevich & Melnick (1981) found the first observational 
confirmation of a correlation between H/3 luminosity and 
line profile width for giant extragalactic HII regions and HII 
galaxies, with residuals that are correlated with the nebular 
metallicity. Subsequent work by Melnick et al. (1987; 1988) 
was devoted to obtain a calibration of this correlation in 
order to make it suitable for distance measurements. 



The distance indicator, defined as M z 



7(o/io 



with O/H the oxygen abundance of the nebular gas, pro- 
vides the predicted luminosity from the relation: 



log L(H fj ) = log M*+P 



(21) 



where the zero-point Po = 29.60 was originally defined from 
a sample of 14 giant extragalactic HII regions (Melnick et al. 
1988) and from which they obtained Ho — 80 ± 5 km sec -1 
Mpc -1 . Obviously, a critical prerequisite for using such seal- 
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ing relations as distance estimators is an accurate calibra- 
tion of their zero-points. Note that a semi-empirical upper 
limit of o = 65 km/s has been proposed for suitable galaxies, 
which can be explained by the requirement that HII galax- 
ies are powered by clusters of coeval starbursts with their 
dynamics dominated by pressure and not rotation. 

The L(Hp) — a relation has been shown to hold also at 
large redshifts (Koo et al. 1996, Pettini et al. 2001, Erb et al. 
2003) . Melnick et al. (2000) showed that HH-like starburst 
galaxies up to z ~ 3, satisfy the L(H/3)-a relation, opening 
the possibility of using the relation to measure cosmological 
parameters. They derived the following distance modulus 
relation for HII galaxies: 

fi = 2.5 log(a 5 /F H 0) - 2.5 \og(0/H) - A H /3 + Z , (22) 

where Fjjp and Anp are the flux and extinction in Hp, re- 
spectively. The originally determined zero-point was Zq = 
—26.18 and the rms distance modulus dispersion was found 
to be ~ 0.52 mag. Although, such an rms uncertainty is 
larger than what is obtained with SNIa, the advantage of 
using HII galaxies is that we can reach a much larger red- 
shift limit (z ~ 4 vs z ~ 1.7). 

Using recent galaxy distance determinations we should 
be able to better determine the zero-point of the distance 
indicator, Zq. To this end we have repeated the original 
analysis of Melnick et al. (1988; 2000), using Cepheid and 
RRLyrae distance determinations and indeed the rms scat- 
ter of the distance indicator relation is reduced by ~ 7% 
while Po = 29.44. This results in a reduction by ~ 0.42 
mags of the zero-point (ie., Zo = —26.60), which provides 
results consistent with Ho = 73 km sec -1 Mpc -1 (Chavez 
et al. 2012 in preparation). 

It should also be mentioned that there are some sys- 
tematic effects than can bias distances obtained with the 
L(Hp) — a relation, in particular differences in the ages of 
the stellar populations, contamination from underlying old 
stellar components, or different extinction laws. To some ex- 
tent these effects can be mitigated by using the equivalent 
widths of the lines to select only very young objects, and the 
use of modern instrumentation that allows a precise control 
of the size, orientation, and location of the spectrograph 
slits. Still the observations remain challenging and require a 
high level of planning and control. 

We are at the process of completing an investigation 
of these effects by using high-resolution spectroscopy of a 
relatively large number of SDSS low-2 HII galaxies with a 
wide range of relevant parameters (Hp equivalent widths 
and luminosities, metal content, and local overdensity) in an 
attempt to understand systematics and to reduce the scatter 
of the distance estimator (Chavez et al 2011 in preparation). 

Most high-z HII galaxies known until recently were 
found in broad-band searches aimed mostly to search for 
Lyman break galaxies, which means that they generally 
have relatively strong continua. Still, a substantial fraction 
present strong emission lines making them ideal for our dis- 
tance estimator (see for example Erb et al. 2006a; 2006b) 
Furthermore, deep slit-less surveys using WFC3 on HST and 
Narrow band filters at SUBARU have revealed substantial 
numbers of HII galaxies with large equivalent widths (i.e. 
strong emission lines and weak continua) at intermediate 
and high redshifts (Yamada et al. 2005; Kakazu, Cowie, & 
Hu 2007; Xia et al. 2011; Atek et al. 2010; Nestor et al. 2011; 



Table 2. Cosmological parameter fits using the Siegel et al. HII- 
galaxies and the newly derived zero-point Zq (eq. 22). The QDE 
equation of state parameter, w, remains completely unconstrained 
by the current analysis. 
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Figure 10. The Hll-galaxy QDE constraints (in the f! m , w 
plane), based on the Siegel et al. sample after excluding two HII 
galaxies showing strong indications for a rotational velocity com- 
ponent. Although the constraints are weak, leaving completely 
unconstrained the value of w, they are consistent at a ~ Iff level 
with the SNIa results (thin red contours). 



Straughn et al., 2011). In all, the present sample has more 
than 400 HII galaxies covering the redshift range 0.5<z<3.7 
with about 100 in the range 3.0<z<3.7 and about 150 at 
z~2. 

Summarizing, the use of HII galaxies to trace the Hub- 
ble relation, as an alternative to the traditionally used SN 
la, is based on the following facts: 

(a) local and high z HII galaxies define a phenomenologi- 
cal relation between Hp luminosity, velocity dispersion, and 
metallicity as traced by O/H that holds out to cosmologi- 
cal distances. Thus, HII galaxies can be used as alternative 
tracers of the Hubble expansion; 

(b) HII galaxies can be readily observed at much larger red- 
shifts than those currently probed by SNIa samples, and; 

(c) it is at such higher redshifts that the differences between 
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the predictions of the different cosmological models appear 
more vividly (see Fig. f). 

A more recent application using 15 starburst galaxies 
with redshift in the range z = 2.17-3.39 has been carried out 
by Siegel et al. (2005) in an attempt to constrain cosmo- 
logical parameters, but the resulting constraints are rather 
weak. They found Q m = 0.2llg'j2 f° r a A-dominated Uni- 
verse. Clearly, the errors are still large, and up to this point 
the results are not competitive with SNIa. 

We have performed our own re-analysis of this data-set, 
following a similar procedure to that applied to the SNIa 
data in the previous sections, but using the newly derived 
value of Zo in eq. (22), allowing for the asymmetric uncer- 
tainties of the HH-galaxy distance moduli and correcting for 
the effects of gravitational lensing, according to section 2.2. 
Furthermore, we have updated the values of the stellar ve- 
locity dispersion and its uncertainty, for some of the galaxies 
in the Siegel et al. sample, according to Erb et al. (2006a). 

The resulting constraints on the (fl m ,w) plane are indi- 
cated in the first two rows of Table 2, while in the last two 
rows (and in Figure 10) we present results after excluding 
two HII galaxies (Q1700-MD103 and SSA22a-MD41) that 
show indications of a significant rotational velocity com- 
ponent (derived from the tilted emission lines; Erb et al. 
2006a), which contaminates the estimate of their velocity 
dispersion. 

Our results show that the derived Q, m values, indepen- 
dent of the value of w, are towards the lower end of the gen- 
erally accepted range, while when excluding the two rotat- 
ing galaxies the fitted Q m parameter moves towards higher 
values, while there is also a decrease of the value of the cor- 
responding reduced \ 2 - 

In any case, the main qualitative result of our HII- 
galaxy based analysis is that although the constraints in the 
(f2 m , w) plane are consistent with those of the Constitution 
SNIa analysis, as can be seen in Fig. 10, the provided cosmo- 
logical parameter uncertainties are significantly larger and 
the degeneracy between fl m and w is even more exacerbated. 
These results clearly indicate that the distance indicator for 
HII galaxies is highly competitive provided: 

• We increase the number of available high- 2 (ztt 2) HII 
galaxies, but it is important to cover also the 0.2 < z < 1 
range (as shown in section 2.3.4). 

• We apply the estimator to a significant sample of bona- 
fide high redshift HII galaxies selected by the strength of 
their emission lines to ensure no contamination by rotation 
and/or underlying old stellar populations. 

• We minimize all possible sources of systematic and ran- 
dom errors. 



3 COSMOLOGICAL PARAMETERS FROM 
THE CLUSTERING OF X-RAY AGN 

The method used to put cosmological constraints, based on 
the clustering of some extragalactic mass-tracer (Matsub- 
ara 2004, Basilakos & Plionis 2009 - hereafter BP09 - and 
references therein), consists in comparing the observed spa- 
tial or angular clustering with that predicted by different 
primordial fluctuations power-spectra, in the latter case us- 
ing also Limber's integral equation (Limber 1953) to invert 



the spatial to angular clustering. By minimizing the differ- 
ences of the observed and predicted correlation function, 
one can constrain the cosmological parameters that enter 
in the power-spectrum determination as well as in Limber's 
inversion. Using the latter we can relate the angular and 
spatial clustering of any extragalactic population under the 
assumption of power-law correlations and the small angle 
approximation. 

We have chosen X-ray selected AGN as a tracer of the 
large-scale structure, in order to perform the previously de- 
scribed analysis, for the following reasons: 

(a) X-ray selected AGN can be detected out to high red- 
shifts (the peak of their ^-distribution is ~ 1) and thus trace 
the distant density fluctuations providing a further anchor 
of the evolution parameter at a redshift other than z ~ 0, 
which most galaxy samples trace to-date. 

(b) AGN selected through their X-ray emission (and not in 
the optical) provide a relatively unbiased census of the AGN 
phenomenon, since obscured AGN, largely missed in optical 
surveys, are included in X-ray surveys. 

(c) Furthermore, determining the clustering at (z) ~ 1 and 
z ~0, one can put better constraints on the cosmic evolution 
of the AGN phenomenon and the evolution of the relation 
between AGN activity and Dark Matter (DM) halo hosts 
(eg. Mo & White 1996, Sheth et al. 2001), and finally also on 
the cosmological parameters and the dark-energy equation 
of state (eg. Basilakos & Plionis 2005; 2006; 2009; 2010). 

3.1 Clustering of X-ray AGN: Biases and 
Systematics 

The earlier ROSAT-based analyses (eg. Boyle & Mo 1993; 
Vikhlinin & Forman 1995; Carrera et al. 1998; Akylas, Geor- 
gantopoulos & Plionis, 2000; Mullis et al. 2004) provided 
conflicting results on the nature and amplitude of high z 
AGN clustering. With the advent of the XMM and Chandra 
X-ray observatories, many groups have attempted to set- 
tle this issue, but in vain. Different surveys have provided 
again a multitude of conflicting results, intensifying the de- 
bate (eg. Yang et al. 2003; Manners et al. 2003; Basilakos 
ct al. 2004; Gilli et al. 2005; Basilakos et al 2005; Yang et 
al. 2006; Puccetti et al. 2006; Miyaji et al. 2007; Gandhi et 
al. 2006; Carrera et al. 2007; Coil et al. 2009; Starikova et 
al. 2010). However, strong indications exist for a flux- limit 
dependent clustering, interpreted as an X-ray luminosity de- 
pendent clustering, which appears to remove most of the 
above inconsistencies (Plionis et al. 2008). Such a luminos- 
ity dependent clustering trend was recently reported also by 
Cappelluti et al. (2010) and Krumpe, Miyaji & Coil (2010). 

Furthermore, there are indications for a quite large high 
z AGN clustering length, reaching values ^ 10 ft -1 Mpc 
at the brightest flux-limits (eg., Basilakos et al 2004; 2005, 
Puccetti et al. 2006, Plionis et al. 2008; Cappelluti et al. 
2010), which, if verified, has important consequences for the 
AGN bias evolution and therefore for the evolution of the 
AGN phenomenon (eg. Miyaji et al. 2007; Basilakos, Plio- 
nis & Ragone-Figueroa 2008 - BPR08 hereafter). An inde- 
pendent test of these results would be to establish that the 
environment of high z AGN is associated with large DM 
haloes, which being massive should be more clustered (work 
in progress). 
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It is also important to understand and overcome the 
shortcomings and problems that one is facing in order to 
reliably and unambiguously determine the clustering prop- 
erties of the X-ray selected AGN. Such a list of problems 
includes the effects of Cosmic Variance, the so-called ampli- 
fication bias, the reliability of the log N — log S distribution 
of the X-ray AGN luminosity function, etc. (see discussion 
in Plionis et al. 2009). 

Recently, Ebrero et al. (2009a) derived the angular cor- 
relation function of the soft (0.5-2 keV) X-ray sources using 
1063 XMM-Newton observations at high galactic latitudes 
(2XMM survey). A full description of the data reduction, 
source detection and flux estimation are presented in Ma- 
teos et al. (2008). Note, that the survey contains ~ 30, 000 
soft-band point sources within an effective area of ~ 125.5 
deg 2 (for f x > 1.4x 10 -15 erg cm~ 2 s _1 ). The large area cov- 
ered and the corresponding large number of X-ray sources 
ensure that the previously mentioned cosmic variance ef- 
fects are minimized. However, further details regarding the 
various biases that should be taken into account (the ampli- 
fication bias and integral constraint), the survey luminosity 
and selection functions as well as issues related to possible 
non-AGN contamination, which are estimated to be 10%, 
can be found in Ebrero et al (2009b). 

3.2 Cosmology from the 2XMM angular 
Clustering 

An optimal approach to unambiguously determine the clus- 
tering pattern of X-ray selected AGN would be to determine 
both the angular and spatial clustering pattern. The reason 
being that various systematic effects or uncertainties enter 
differently in the two types of analyses. On the one side, us- 
ing the angular two-point correlation function, w(9), and its 
Limber inversion, one bypasses the effects of redshift-space 
distortions and uncertainties related to possible misidenti- 
fication of the optical counter-parts of X-ray sources. On 
the other side, using spectroscopic or accurate photometric 
redshifts to measure the spatial, £(r), or projected, w p (6), 
2-point correlation function one by-passes the inherent ne- 
cessity, in Limber's inversion of w(6), of assuming a source 
redshift-selection function (for the determination of which 
one uses the integrated X-ray source luminosity function, 
different models of which exist). 

The basic integral equation relating the angular and 
spatial correlation functions is: 

= {~) 2 E(z)*z£«r,zydu, (23) 

where dN/dz is the source redshift distribution, estimated 
by integrating the appropriate source luminosity function 
(in our case that of Ebrero et al. 2009b), folding in also 
the area curve of the survey. Note that to derive the spatial 
correlation length from eq. (23), it is necessary to model the 
spatial correlation function as a power law, assume the small 
angle approximation as well as a cosmological background 
model. The latter is provided by the function E(z) (eq. 4), 
which for a fiat background and the QDM equation of state, 
it takes the form: 

E(z) = [fl m (l + zf + (1 - fi m )(l + z)'^ 1+w) ] 1/2 . (24) 
The AGN spatial correlation function can be written as: 



ar,z) = (l + z)-^b 2 (z)^M(r) , (25) 

where b{z) is the evolution of the linear bias factor (eg. Mo 
& White 1997; Matarrese et al. 1997; Sheth & Tormen 1999; 
Basilakos & Plionis 2001; 2003, Basilakos et al. 2008; Tinker 
et al. 2010; Ma et al. 2011) e is a parameter related to the 
model of AGN clustering evolution (eg. de Zotti et al. 1990) 
and £dm(j") is the corresponding correlation function of the 
underlying dark matter distribution, given by the Fourier 
transform of the spatial power spectrum P(k) of the matter 
fluctuations, linearly extrapolated to the present epoch: 

The CDM power spectrum is given by: P(k) = 
P k n T 2 (Q. m ,k), with T(Q m ,k) the CDM transfer function 
(Bardeen et al. 1986; Sugiyama 1995), n ~ 0.967 and a 
baryonic density of Q^h 2 = 0.02249, following the 7-year 
WMAP results (Komatsu et al. 2011). The normalization 
of the power-spectrum, Po, can be parametrized by the rms 
mass fluctuations on Rs = 8ft _1 Mpc scales (as), according 
to: 

Po = 2n 2 a 2 /9(fl m ,R s ) , (27) 
with 

V(n m ,Rs)= k n+2 T 2 (Q m ,k)W 2 (kR 8 )dk (28) 
Jo 

and W(kR 8 ) = 3(sinfci? 8 - kR$coskR$)/(kR 8 ) 3 . 

Evidently, the essential parameters needed to charac- 
terize any QDE cosmological model are: Q m ,w,a 8 and Ho- 
Regarding the Hubble constant we will use the WMAP7 re- 
sults (Komatsu et al. 2011), which practically coincide with 
those of the HST key project (Freedman et al. 2001), ie., 
h — Ho/100 = 0.704, while regarding the as normalization 
of the CDM power spectrum we will use the extrema of the 
range provided by the recent analysis of SDSS LRGs for a 
range of dark energy equations of state (a b € [0.78,0.81]; 
Sanchez et al. 2009). Note that the upper limit of the above 
range corresponds to the WMAP7 ACDM value (Komatsu 
et al. 2011). 

Furthermore, to estimate the predicted QDE model cor- 
relation function of the underlying mass, £dm(>", z), in order 
to compare it with the observed AGN clustering, it is nec- 
essary to deal with the following three issues: 

(i) Clustering Evolution Model: As discussed earlier 
(see eq. 25), in order to estimate the expected clustering 
of any mass tracer it is important to assume a clustering 
evolution model (eg. de Zotti et al. 1990), which is encap- 
sulated in the value of the parameter e. A value e = —1.2 
corresponds to a constant in comoving coordinates cluster- 
ing model, while a value e = —3, to a constant in physical 
coordinates. According to Kiindic (1997) and Basilakos & 
Plionis (2005; 2006) we will use the former value of e (al- 
though we have also tested the effects of using e = —3). 

(ii) Bias Evolution Model: We need to calibrate the 
parameters of the bias evolution model to each cosmologi- 
cal model. Although a large number of bias evolution mod- 
els have been proposed in the literature (see Papageorgiou, 
Plionis & Basilakos 2011, in preparation for a comparison 
of different models) , we use here the approach of Basilakos 
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& Plionis (2001; 2003) which was extended to QDE cos- 
mological models in BPR08. This model is based on linear 
perturbation theory and the Friedmann-Lemaitre solutions 
of the cosmological field equations, and includes also the 
effects of interactions and merging of the mass tracers. Its 
analytical form has been derived for the QDE cosmological 
models, and its generalization to the CPL and alternative 
gravity cosmological models is under-way (Basilakos, et al. 
in preparation). Considering that each X-ray AGN is hosted 
by a dark matter halo of mass Mh, we can analytically pre- 
dict its bias evolution behaviour within the QDE models. 
Conversely fitting the model to observations, we can deter- 
mine the mass of the DM halo within which AGN live (for 
more details see Basilakos & Plionis 2010). 

For the case of a spatially flat cosmological model, our 
bias evolution model has the following form: 

b(M h ,z) = C 1 (M h )E(z) + C 2 (M h )E(z)I(z) + y p (z) + l ,(29) 

where y P (z) determines the rate of halo merging^. However 
it is important only for 3 and therefore we neglect it here 
(see BPR08). Furthermore, we have 



I{z) = 



( 1 + - 



\E(a 



dx . 



(30) 



while the constants Ci and Ci have been fitted using a 
ACDM simulation with « m = 0.3 and a A = 0.9 (BPR08), 
and have been found to follow the form: 

CAM h ) * ai , 2 [j^j^ ) (3D 

with Pi = 0.34, /?2 = 0.32, while ot\ and a.2 have been found 
to be cosmological model dependent, with values given by 
(Papageorgiou, Plionis & Basilakos, 2011 in preparation): 

'0.! 



Ki 



exp [K 3 (^m - 0.3)] 



with Ki — 3.44, «2 — 2/5 and K3 ~ 4/5, and 



Ct2 



-0.36 



I 0.3/ 



3/2 



(32) 



(33) 



(iii) Non-linear Power Spectrum: Since the correla- 
tion function on small angular scales is within the expected 
non-linear regime, we should include in our model power- 
spectrum the non-linear contributions. To this end we use 
the corresponding fitting formula introduced by Peacock 
& Dodds (1996), for the ACDM model (see also Smith et 
al. 2003; Widrow et al. 2009). There is one relatively free 
parameter in their formulation, which is the slope of the 
power spectrum at the relevant scales, since the CDM power- 
spectrum curves slowly and thus it varies as a function of 
scale: n c g — din P(fc)/dln k. On the scales of interest the 
value is n e a — —2, but we have decided to actually derive, 
and then fix, the n e fr value from the data analysis itself. 
Using the minimization procedure discussed in section 3.3 
we compare the observed 2XMM AGN correlation function 
with that provided by the WMAP7 ACDM model (ie., fixing 
Q m = 0.272, w= — 1, and erg = 0.811) fitting the remaining 
two free parameters (ie., Mh and n e ff). The corresponding 



S Note that the bias factor at the present time is given by: 
b(M h ,0) = Ci(M h )+C 2 (M h )I(0) + 1. 




M h /10 13 h- 1 M s 

Figure 11. Main Panel: The 1, 2 and 3<r likelihood contours in 
the M h ,n c {{ parameter space for the ACDM (WMAP7) model. 
Inner Panel: The Ebrero et al. (2009a) 2XMM angular correlation 
function and the best fit ACDM model (continuous line). 



solution space can be seen in Fig. 12, while the best fit 
parameter values are: Mh ~ 6.5(±2.1) x 10 12 ft -1 Mq and 
n c g ~ —2.02^004. In the inset of Fig. 11 we also plot the 
2XMM angular correlation function together with the best 
fit ACDM model (continuous line). 

As a consistency check we have verified that when fixing 
the non-linear slope of P(k) to the above fitted value and 
leave as free parameters Mh and as , we recover the WMPA7 
as value and exactly the same Mh value, as above. The 
derived value of Mh is slightly larger than that provided 
by Ebrero et al. (2009a) using the Sheth et al. (2001) bias 
evolution model (ie., ~ 5 x 10 12 /i _1 Mq). We have tested 
also the case of a clustering evolution model with e — —3, 
in which case the derived value of the halo mass is ~ 2 x 
10 10 /i -1 Mq, a value significantly below any reasonable value 
that has been proposed or derived in the literature. We will 
therefore use e = —1.2 throughout the rest of the paper. 

3.3 Fitting Models to the 2XMM Clustering Data 

In order to constrain the cosmological parameters we use 
again the standard x 2 likelihood procedure and compare the 
measured XMM soft-band source angular correlation func- 
tion (Ebrero et al. 2009a) with the predictions of different 
spatially flat cosmological models. The corresponding likeli- 
hood estimator is defined as: £agn(p) oc exp[— Xagn(p)/2] 
with: 



Xagn(p) 



E 

»=r 



[tuth(0i,p) - w ohs (8i)] 2 



(34) 



where p = (fi m , w, Mh), c< is the uncertainty of the ob- 
served angular correlation function and ae i corresponds to 
the width of the angular separation bins. 

We sample the various parameters as follows: the matter 
density f2 m £ [0.1, 0.4] in steps of 0.002; the equation of state 
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parameter w£ [—1.4, —0.6] in steps of 0.005 and the parent 
dark matter halo M h /W 13 h~ 1 M e , G [0.1, 3] in steps of 0.01. 
In this likelihood analysis we use as priors a fiat universe, 
and the previously mentioned values of h, as and 

The results of the minimization procedure for the case 
of (j 8 = 0.81 are: ft m = 0.301 ± 0.008, w= -0.990 ± 0.058 
and M h = 3.1(±1.1) x 10 12 h' 1 M , with a X 2 = 39.41 for 
10 degrees of freedom. The large value of x 2 /dof is due to 
the sinusoidal modulation of the 2XMM w(6), which could 
be due to systematic effects possibly related to the size of 
the XMM fields (see discussion in BP09). If we use the 2a 
w(8) uncertainty in the denominator of the x 2 function of 
eq. (34), then the \ 2 drops to ~ 9.86 (and the uncertainties 
of the fitted parameters increase to roughly twice the values 
indicated previously). 

These results slightly differ with respect to the similar 
analysis of BP09, due to a number of improvements that we 
have currently included, apart from the fact that we have 
also used the WMAP7 cosmological parameters (ie., Qb, n 
and h). The two main improvements have to do with (a) 
the bias evolution model, in which we have now taken into 
account the dependence of the parameter ai (see eq.32) on 
Q. m and as (which has mostly affected the derived value of 
Mh, reducing it significantly), and (b) the non- linear power 
spectrum corrections, for which we have used the Peacock 
& Dodds (1996) ACDM fitting formula. 

Nevertheless our current procedure can and will be im- 
proved in the future in a number of ways: 

• We will eventually use the clustering of X-ray selected 
AGN from a large contiguous X-ray survey, a fact which will 
solve, we believe, the quasi-sinusoidal small amplitude mod- 
ulation of the 2XMM w(9) (see discussion in BP09). Such a 
future survey is the XMM-Newton Very Large Programme 
(XXL), which was recently granted time to map two extra- 
galactic regions of 25 deg 2 , at a depth of ~ 5 x 10~ 15 erg 
cm" 2 s _1 (Pierre et al. 2011). 

• We will investigate more accurate non-linear power- 
spectrum corrections of w(8) (eg., Widrow et al. 2009), and 

• We will ultimately test a large range of Dark Energy 
equations of state, to include CPL and alternative gravity 
[f(R)] models. To this end we will use the recent generaliza- 
tion of the BPR08 bias evolution model of Basilakos, Plionis 
& Pouri (2011). 



4 JOINT HUBBLE- RELATION, AGN 
CLUSTERING AND BAO ANALYSIS 

Here we will perform an example of the joint analysis be- 
tween the previously discussed results from 2XMM clus- 
tering and the Hubble relation, which is the basic aim of 
our overall project. For the current exercise we will use the 
Constitution SNIa Hubble relation since we are still working 
on the development of the HH-galaxy based methodology. 
For completeness we will also use the recent results of the 
baryonic acoustic oscillation technique (BAO). We remind 
the reader that BAOs are produced by pressure (acoustic) 
waves in the photon-baryon plasma in the early universe 
generated by dark matter overdensities. At the recombina- 
tion era (z ~ 1100), photons decouple from baryons and free 
stream while the pressure wave stalls. Its frozen scale, which 
constitutes a standard ruler, is equal to the sound horizon 



length, r s ~ 100 h' 1 Mpc (e.g. Eisenstein & Hu 1998). This 
appears as a small, ~ 10% excess in the galaxy, cluster, 
or AGN power spectrum (and its Fourier transform, the 2- 
point correlation function) at a scale corresponding to r s . 
First evidence of such an excess have been reported in the 
clustering of the SDSS luminous red galaxies (LRGs) (see 
Eisenstein et al. 2005, Padmanabhan et al. 2007; Percival 
et al. 2010). In this work we use the latest measurement of 
Percival et al. (2010): 

r s (z d )/D v (z*) = 0.1390 ± 0.0037 , 

(see also Kazin et al. 2010a, 2010b). Note that r s (z d ) is 
the comoving sound horizon size at the baryon drag epoch 
Zd, which is given by the fitting formula of Eisenstein & 
Hu (1998), Dv(z+) is the effective distance measure and 
z* = 0.275. Of course, the quantities r s ,Dv can be defined 
analytically, and are given by: 

_c_ [ a * do_ 

v3. 



,{Zd) 



(35) 



a 2 H (a) sj\ + (3ft(,/4fi 7 )a 

where a d = (1 + Zd)^ 1 and fi 7 /i 2 ~ 2.47 x 10~ 5 the energy 
density of photons. In this context, the effective distance is 
(Eisenstein et al. 2005): 

1 1/3 



D v (z) 



(l + z) 2 D 2 A (z) 



H(2 



(36) 



where Da(z) is the angular diameter distance. Therefore 
the corresponding Xbao function is simply written as: 

r ^^-(p)- 0.1390] 2 



2 / n lD v ( Zi ,) 

Xbao(P) = — 



0.0037 2 (37) 
where p is the vector containing the cosmological parameters 
that we want to fit for. In this case p = (f2 m ,w). 

We therefore perform a joint likelihood analysis, assum- 
ing that any two pairs of cosmological data sets are indepen- 
dent (which indeed they are) and thus the joint likelihood 
can be written as the product of the two individual ones. 
The results based on the joint analysis of the different pairs 
of cosmological data are shown in Fig. 12 and quantified in 
Table 3. It is evident that the addition of the XMM cluster- 
ing analysis provides significantly more stringent constraints 
than, for example, the joint SNIa and BAO results. The re- 
duction parameter, ie., the ratio of the Figure of Merit of the 
joint XMM-SNIa analysis to that of the Constitution SNIa 
analysis (see definition in section 2.3.4) shows that the 2a 
range of the (fi m ,w) solution space is reduced by a factor 
of ~5/2 with respect to that of the BAO-SNIa analysis (see 
lower right panel of Fig. 12). 

The necessity, however, to impose constraints on a more 
general, time-evolving, dark-energy equation of state (eq. 2), 
implies that there is ample space for improving the current 
analysis and indeed our aim is to develop further this project 
by (a) using a new Hubble relation analysis, based on high 
z HII galaxies, as detailed in this paper, and (b) by general- 
izing the BPR08 bias evolution model for any DE equation 
of state (CPL and alternative gravity models). 



5 CONCLUSION 

We have investigated the question of which is the most ef- 
ficient strategy to tighten the cosmological constraints pro- 
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Figure 12. QDE model parameter constraints (ie., in the f2 m ,w plane) provided by the joint likelihood analysis of pairs of cosmological 
data. The red shaded contours are the joint-likelihood contours of the indicated pairs of data. For clarity we show only contours 
corresponding to the 1 and 3 a confidence levels. The 2XMM results shown correspond to those based on a power spectrum normalization 
of a% = 0.81. Note also that we have followed the conservative approach of using 2a w(9) uncertainties in the minimization process, to 
allow for the small amplitude sinusoidal w(9) modulation (see main text and BP09). Upper Left Panel: Constitution SNIa Hubble relation 
(blue contours) and 2XMM AGN clustering (black contours). Upper Right Panel: Constitution SNIa Hubble relation (blue contours) and 
LRGs BAO (green contours). Lower Left Panel: 2XMM AGN clustering (black contours) and LRGs BAO (green contours). Lower Right 
Panel: The joint likelihood contours of the SNIa-2XMM (red contours) and of the SNIa-BAO (blue contours) pairs. 



vided by fitting the Hubble relation. Using extensive Monte- 
Carlo simulations we have verified that by using only a 
few high z tracers (in the range 2^$ z^, 3.5), even with a 
relatively large distance modulus uncertainty, we can re- 
duce significantly the present cosmological parameter solu- 
tion space. We have taken into account the effects of lensing 
magnification/de- magnification, which not only increases 
the distance modulus uncertainty but it also shifts system- 
atically the mean distance modulus of individual sources. 
Although the effects can be severe for an individual source, 



they can be statistically treated and they are significantly 
reduced the denser the source sampling is in redshift space. 
Applying our lensing magnification correction to the Con- 
stitution SNIa set (Hicken et al.), we find that the fitted cos- 
mological parameters are not significantlly affected by such 
effects, due to the fact that the SNIa sample traces relatively 
small redshifts (z^S 1). 

Based on a figure of merit analysis we have provided 
a simple procedure to estimate the necessary number of 
2;$ z^ 3.5 tracers needed to reduce the cosmological solu- 
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Table 3. The best fit values of the cosmological parameters based on the joint likelihood analysis of the indicated cosmological data. 
For the case of the 2XMM clustering analysis we followed the conservative approach of using 2<r w(8) uncertainties. The uncertainty of 
each fitted cosmological parameter has been estimated after marginalizing over the other parameter (ie., by fixing one parameter at its 
best value and allowing the other to vary, providing as its uncertainty the range for which A% 2 < 2.3). The last column indicates the 
reduction parameter (related to the Figure of Merit) as defined in section 2.3.4. 



Joint data 




w 


x 2 


dof 


S 


SNIa+BAO 


0.296 ±0.021 


-1.027 ±0.053 


439.711 


365 


2.32 


XMM a +SNIa 
XMM"+BAO 


0.310 ±0.012 
0.302 ±0.014 


-1.064±°;°il 
-0.995±8;°g| 


449.591 
9.860 


374 
11 


4.49 
1.74 


XMM 6 +SNIa 
XMM<>+BAO 


0.318 ±0.013 
0.306 ±0.015 


-1.085±g;8S 
-0.973±°;°H 


449.644 
9.861 


374 
11 


4.35 
1.70 



a 2XMM results based on a P(k) normalization of erg = 0.811, 
b based on cr 8 = 0.78. 



tion space, presently provided by the Constitution set, by 
a desired factor of our choice and for any level of rms dis- 
tance modulus uncertainty. This analysis has shown that in 
order to significantly reduce the cosmological parameter so- 
lution space, it is more efficient to increase the number of 
high- z tracers than to reduce their individual uncertainties. 
A re-analysis of the cosmological constraints provided by a 
small sample of high-z (2.2< z< 3.4) HII galaxies, previ- 
ously analysed by Siegel et al. (2005), but now using a novel 
determination of the zero-point of the relevant distance scal- 
ing relation, provides consistent cosmological results with 
those of SNIa, although the dark energy EoS parameter re- 
mains unconstrained at present. 

Finally, using the clustering of X-ray selected AGN we 
provide the framework that will be used, joining their cosmo- 
logical likelihood with that of the Hubble relation analysis, 
to put stringent dark energy equation of state constraints. 
An example of such a joint analysis, using the 2XMM clus- 
tering and the Constitution SNIa Hubble relation, and under 
the priors of a flat universe, h = 0.704 and a$ — 0.81 or 0.78, 
provide significantly more stringent QDE model constraints, 
as indicated by the fact that the Figure of Merit increases 
by a factor ~ 2, with respect to that of the joint SNIa- 
BAO analysis. The QDE cosmological parameters provided 
by the 2XMM-SNIa joint analysis are: f2 m = 0.31 ±0.01 and 
w= —1.06 ± 0.05, with the uncertainties being estimated 
after marginalizing one parameter over the other. 
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